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Abstract 

We consider random walk on a mildly random environment on finite transitive d- 
regular graphs of increasing girth. After scaling and centering, the analytic spectrum 
of the transition matrix converges in distribution to a Gaussian noise. An interesting 
phenomenon occurs at d = 2: as the limit graph changes from a regular tree to the 
integers, the noise becomes localized. 
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1 Introduction 

Localization phenomena for eigenvalues of random media have received considerable atten- 
tion lately. Among the several results, we point out a specific one. If one perturbs the 
Laplacian of the nearest-neighbor graph on the integers by a very small i.i.d. potential, its 



spectrum becomes pure point immediately. In cont rast, t 
regular tree preserves the continuous spectrum, see 



Klein 



lesame change on a higher degree 



The results of this paper point to an analogous, but distinct phenomenon. While the 
above localization phenomena are due to large-scale behavior of the eigenvectors, we find a 
similar dichotomy for the local behavior. 

For this, we consider sequences of vertex-transitive finite graphs G„ with degree d > 2 
and increasing girth (the length of the shortest simple cycle) ; these converge locally to either 
the integers {d = 2) or a regular tree of higher degree [d > 3). On these graphs, we 
consider a perturbation of the Laplacian, given by the w ell-studie d rand om walk on random 



environment model with small noise (see, for example IZeitounil ( l2004j ). The noise in the 
Laplacian creates a Gaussian noise in the spectrum. Our main discovery is that in the limit 
this noise is local for the integers and has long-range correlations for higher-degree trees. 

We first consider more general graphs G. We will have the following standing assump- 
tions. 

Assumptions. G is a vertex-transitive graph of finite degree. Let M denote the tran- 
sition probability matrix for nearest-neighbor simple random walk on G. For the random 
environment on G one randomly modifies transition probabilities along edges. Its transition 
matrix is defined as M + eB, where i? is a random matrix satisfying the following. 

The random variables By_^y have mean zero and variance 1. Further, the B^^^^s corre- 
sponding to different u are independent, and for some constant Ci and all vertices u,v E G 
we have 

\Bu,v\ < CiMu,v, (1) 

Y. Bu,,. = 0. (2) 

geStab(u) 

Here Stab('u) is the stabilizer of u in the automorphism group of G. Condition (^ means 
that the sum of the bias of the random environment over symmetric directions is zero, which 
is a bit stronger than just requiring M + eB to be stochastic. Condition ([1]) implies that 
Bu^y vanishes when uv is not an edge and that for small e the entries of i? + eM are in [0, 1]. 
These assumptions will be in effect for the rest of the paper. 

For a finite G, let /i^ denote the empirical probability measure of the eigenvalues of 
M + eB. Then as e ^ we will show an expansion 

where /i" is a random functionals. We define the second difference quotient 

e^ 



where the extra factor upfront makes the scahng consistent as G changes. This is the centered 
and scaled empirical eigenvalue measure of M + eB. Let X denote the space of complex 
functions analytic in a neighborhood of the closed unit disk. We call / real if it maps M to 
M. In Proposition IH we will show that for / G A* the limit 

TG(/) = lim I fdm, (3) 

exists. Our first theorem identifies the covariance structure of Tclf) as / ranges over real 
functions in X. 

For a possibly infinite vertex-transitive graph G and complex |A| < 1, let 

px{x) = 1/(1 - Xx), Gx = Px{M) = {I- \M)-\ 

the Green's function corresponding to G. Let o be a marked vertex of G and define 

Hg{px,P^) := ^dl^ ( AVI]E[(5^a)o.(56^;.)o.]E[(56;a).o(56^m).o]] ) . (4) 

We will show in Lemma [7] and Proposition [8] that Hg{p\,p^) is well-defined and extends 
uniquely to a bilinear form Hq on X. 

Theorem 1. For any finite vertex-transitive graph G and f,g^Xwe have 

ETcif) = 0, E[TGif)TGig)] = Hcif^g). 

We now consider sequences of transitive graphs (?„ -^ G locally, which means that for 
every r the r-neighborhood of a fixed vertex of Gn eventually agrees with that in G. 

Theorem 2 (Normality). Let Gn —>■ G, an infinite graph, and let Gn, G satisfy the assump- 
tions above with uniform constant Ci. Then as n —>■ oo, jointly for all real f ^ X the random 
variables Tc^if) converge weakly to mean zero normal random variables T{f). Moreover, 
for all f,gEXwe have 

E [TgM)TgM] ^ E [T{f)T{g)] = Hdf^g). 

Finally, we consider the case when Gn -^ ^d, the d-regular tree. Since T^^ is bipartite, we 
expect that the random functionals me have a symmetric limit. Therefore it is natural to 
consider T„ for even functions of the form f o s where s{x) = x"^ and f & X. Note that for 
any finite graph G 

TG{fos)=fGif), 



where Tdf) is defined analogously to Tdf) but in terms of /i^, the empirical eigenvalue 
measure of (M + eBY. 



It is known (see 



Linial and Widgersonl . Theorem 13.4) that yU„ corresponding to Gn con- 



verges weakly to the measure with density 

2d'^ ^xi^p — x) 



aAx) 



, Q<x<pi (5) 

TT 1 — X 

where p = pd '■= 4:{d — l)/d'^, the squared spectral radius of the walk on the infinite (i-regular 
tree T^. Our next theorem gives an explicit expression for the limiting covariance structure. 

Theorem 3 (Tree limits). Assume that Gn — ^ T^, the d-regular tree. Then jointly for all 
real f E X we have the convergence in distribution TcXf) ~^ T{f), where the T{f) are 
jointly normal and have mean 0. Moreover, 

rp rP 



nnf)T{g)] 



JO 



f'{x)g'{y)pd{x,y)dxdy, 



where for d > 2 the kernel (3d is given by 

2rf^ {d-2)K{x)K{y) 



Pd{x,y) 



with the semicircle function 
k{x) :- 



7r2 16{2d -3){x- yf + {d - 2yA{d, x, y) 



(6) 



(7) 



2d\/x{p — 



X] 



and 



A{d,x,y) := pn{^yy + A{d + 3)ix-yf + p\d-2f 



For d = 2, the covariance is given by 



32 



7r 



f'{x)g'{x)\/x{l — x)dx, 

which corresponds to (32{x,y) = —^Jx{l — x) 6x{y)- 

The kernel function jSd becomes singular as rf | 2. This can be seen in the figure on 
the first page showing graphs j3d '■ [0,p]^ —>■ M for d = 4,3,2.1 from above. In words, 
in the d = 2 case the centered and scaled empirical measure converges to the so-called 
i/i-noise with a certain density. In contrast with the ifi/2-noise that arises in the limit 
of the Gaussian Ortho gonal Ensemble and Haar u nitary random matrix n iodels (among 



J2OO1I ) 



others; see, for example iDiaconis and EvansI (I2OOII ). lAnderson and Zeitounil ( 20061 ) and the 
references therein), this noise is local, as the ^-funct ion in the covariance for mula shows. Hi 
noise appears typically for complex eigenvalues, see iRider and Virag) (120071 ). Formula (EI) is 



written in a form to show how the kernel becomes singular (localized) as d I 2. Recall that 
'K^^e/{x'^ + e'^) converges as e — > to the delta function at x = 0. 

The presence of /', g' in the formula for the covariance is a priori expected. When 
the matrix is randomly perturbed, the eigenvalues move in random directions. For every 
eigenvalue A, the function /(A) changes by a small amount proportional to /' times the 
change in A. However, the picture here is more complicated because the perturbation is of 
order e^. The kernel f3 reflects the correlation between the random change in eigenvalues as 
well as their density. 

Theorem [3] motivates the following question: 

Question 1. Which natural sequences of graphs give rise to localized noise? 

Our strategy is as follows. To prove normality, we use a central limit theorem for depen- 
dent variables based on Stein's method, rather than the usual method of moment computa- 
tions. For computing the covariance, we consider traces, and so we will have to count certain 
paths. We use the Green function to jointly treat paths of different lengths and to avoid 
complicated computations with orthogonal polynomials that normally arise in this context. 
Theorems [H [2] and [3] are proved in Sections HI [5] and [7| respectively. 

The finite graphs that we study arise naturally. The simplest way to pick an interesting 
transitive graph i s to consider a C ayley graph of a sufficiently complicated finite group. As 
it is discussed in iGamburd et al.l . such Cayley graphs will typically have girth tending to 
infinity; for example, random Cayley graphs of simple groups of increasing order will have 
this property. Of course, for d = 2, the only examples are graphs consisting of cycles. 

Rando m perturbations of eigen v alues have been extensively studied in the literature. 



Reference 



vom Scheldt and Purkert 



(119831 ) studies properties of the eigenvalues of a random 



perturbation of a fixed matrix without the restriction for stochasticity or positivity of the 



matrices involved (see Sections 1.3 and 2.2). Another possib le approac 



1 to our problem is 



via perturbation expansions for eigenvalues (see for example iDeiJ ( 1l99ll ) Sections 6.3, 6.4), 
but this requires more control of the eigenvectors. 

2 The £ ^ limit 

We first show that for finite graphs G the limit ([3]) defining Tdf) exists and identify it. We 
will rely on the Assumptions about the perturbation matrix B, in particular, ([2]). In fact, 
throughout this paper, we only need a simple consequence of ([2]). Namely, for all A; > 1 and 



V E G, we have 

(5M'^),,„ = 0. (8) 

Indeed, for any g G Stab(f ) this equals 

weG weG weG 

Averaging each term over all g G Stab(w) and using ([2]) we get (jS]). 

Proposition 4 (The limit as e — > 0). For any f = Yl S-^"' ^ '^' ^^^ limit Tdf) exists, and 



^ oo 

TG(/) = -Tr(a,,/(M + £5)|,=o) =l]a,TG(z^'). (9) 

j=0 

Moreover, for integers j > we have 



ki+k2+2=j 



Proof. We use the fact (see, for example. Theorem 6.2.8 in iHorn and JohnsonI (Il99ll )) that 
for / an analytic function with radius of convergence r > 0, and for A in the set Air of 
matrices of some fixed dimension and with spectral radius less than r, the power series f{A) 
is absolutely convergent and is analytic as a function of (the entries of) A. Thus Tr f{M+eB) 
is an analytic function of all of its (at most) 21(71^ + 1 variables near M, B, and e = 0. We 
are free to rearrange its absolutely convergent multiple power series expansion in any way 
we like. So the first derivative with respect to e at is the coefficient of the e terms; this is 
a multiple of 

J2 Tr [M'^'BM''^) = J2 Tk{.BM^'+^'') 

ki+k2=j — l ki+k2=j — l 

and each of the summands is zero by (j8]) . We have 

\G\^I^ j fdm, = e-\i:if{M + eB)-i:if{M))-^]^deeTif{M + eB)l=, 

when £ ^ 0, since the first derivative vanishes at e = 0. Now the right hand side equals the 
coefficient of e^ in the expansion. This gives ([9]). 

Since the trace of a product does not change when we cyclically permute the factors, we 
see that 

de Tr((M + eBY) = Tr(9e(M + eBY) = j Tr(5(M + eBy-^). 

Taking another derivative, we get the second claim of the proposition. ■ 

6 



3 Properties of the bilinear form Hq 

The goal of this section is to define the bihnear form Hq introduced in (j4]) and estabhsh its 
continuity properties. 

Recall that X denotes the set of power series centered at zero, and with radius of con- 
vergence more than 1. For a power series /, let [x^]f{x) denote the coefficient of x'^ in the 
expansion of /. Given a finite or infinite graph G (with B and M) as in the introduction, 
we fix a vertex o of G, and re-define the bilinear form on X as 

HgU.9) ■■ = ^ E^^"^^- Wix)[x']9ix) (10) 



2 



a.j[x'-']f'{x)[x^-']g'{x) 



where 



a 



fJ2^ 



11) 



,, .-E Yl yi{G,o,v)Y,{G,o,v 
veG\{o} 

and for any vertices v, w of G, we define 

Y^{G,v,w):= Yl {BM^%^^{BM- 

jl+J2+2=j 

The last sum is always finite and symmetric in v,w. In Proposition [8] we will show that 
this definition agrees with (jlj) for the functions p\. Note that the definition of Hq does not 
depend on o if G is vertex-transitive. For any graph G, the bilinear form Hcifig) clearly 
makes sense for polynomials /, g. To go beyond polynomials, we need to show that the 
infinite sum in (ITOj) is well-defined. For f = Y1 '^kZ^ G X, set 

+ 00 

k=l 

which defines a norm on X. (To be precise, it is a norm on X modulo the constants; however, 
we note that for the purposes of this paper the constant terms for functions in X do not 
matter, so we may thing of X itself as the space of functions modulo the constants.) Note 
that on X this norm is finite. 

Remark 5. On the subspace of X consisting of power series with radius of convergence 

strictly larger that a fixed r > 1, the norm || • ||* is dominated by the supremum norm on 

Dr := {z : \z\ < r}. Indeed, for f{z) = J2k=o'^kz'' in that subspace, Cauchy's inequalities 

give \ak\ <r-'=||/||^^. So that 

r(l + r) 



< 



1)= 



oo 



As a direct consequence, we get 
Lemma 6. Polynomials are \\ ■ \\* - dense in X . 

We are now ready to show that Hq is well-defined and continuous. 
Lemma 7. For f,g^X and any G, the sum giving Hcifig) is absolutely convergent, and 

\HGU,9)\<c{\\fUg\U. 
Proof. Note that since |-B„^| < ciM„^, we have 

Y,\yAG,vM\<c\Y. E Mi]:;'Mi-^' = ci{j-i){M^u<cij{M%„. (12) 

w&G w&Gji+J2+2=j 

In particular, the same bound holds for each term \Yj{G,v,w)\. Therefore 

^ J2 Y,iG,o,v)Y,iG,o,v)< J2 K,.Ml,v< E K^M^,, = M^^^' < 1. 

^ veG\{o} veG\{o} veG\{o},w€G 

The claim now follows by summing over alH, j. ■ 

Proposition 8. Hq satisfies (j4]). Moreover, (jlj) uniquely defines Hq as a \\ ■ \\^,- continuous 
bilinear form. 

Proof. The absolute convergence of the series defining He-, implies 

Using the definition of aij, we write the sum above as 

veG\{o} 

k,l,k',£'>0 

= aV E E{Bgx)^{Bgx)UBg^)oviBg^u 

v£G\{o} 

= aV E ^[{Bgx)ov{Bg^UB[{BgxU{BgM. 

vGG\{o} 

This proves the first part. For the second, it suffices to show that the linear span of {px : 
|A| < 1} is II • 11^,-dense in X; this will be done in the next lemma. ■ 

Lemma 9. The vector space generated by {px '■ |A| < 1} is \\ ■ ||* - dense in X . 



Proof. Because of the density of polynomials in X and Remark [5], it suffices to prove that 
any polynomial P can be approximated in the supremum norm on a disk Dr with r > 1 by 
elements of the linear span {px : |A| < 1}. 

Pick 1 < r < ri, and let C.^^ : {2; : |2;| = ri}. Then 

^(0 



P z) 



27ri Jcr, C- z 



1 f P(C) 



Call Rn{z) the Riemann sum corresponding to an equipartition of the circle with n pieces. 
The sequence of functions Rn is equicontinuous on {z : \z\ < r} and it converges pointwise 
to P{z), thus the convergence is uniform on that set. So that for a given 5 > 0, there is a 
finite linear combination 

As{z) := cipi/(^^{z) + 02^1/^2(2) + . . . + CkPi/Q^{z) 

such that 

\f{z)-As{z)\<5 (13) 

for all 2; G C with \z\ < r. This completes the proof of the lemma. ■ 

Finally, we check that HG{f,g) is continuous in G as well. Recall the definition of local 
convergence of graphs from the introduction. For the following lemma, we use, as usual, the 
assumptions from the introduction for each Gn and G. 

Lemma 10. If Gn — > G locally, then for f,g & X we have HG„{f,g) -^ Hc^f^g)- 

Proof. Note that as n — >■ cxd, the neighborhood of radius z + j of o in G„ stabilizes to look 
like the same neighborhood in the limit graph G. 

For polynomials f,g, and for all large n, we have HcSf.g) = Hc^f^g)- This is because 
Hcnifig) only depends on a neighborhood of the root of radius given by the maximal 
degree of / and g. 

Now we have that the sequence of functions Hc^if, g) is equicontinuous on X^ (by Lemma 
[7]) and converges on a dense set. Thus, by Lemma (TT^ they converge on the entire set to a 
continuous limit. Hg{-, ■) is continuous (Lemma [7]), and this finishes the proof. ■ 

4 The covariance structure 

After establishing some properties of the bilinear form Hq, we show that for finite graphs G 
it gives the covariance structure of Td^f). 

9 



Lemma 11. For any finite vertex-transitive graph G and complex polynomials f,g, we have 

E(Tg(/)) = 0, 
E{TGif)TG{g)) = HG{f,g). 

Proof. We will show this for monomials. The extension to polynomials is straightforward 
from bilinearity. First, we have 

Tg(z^) = -^=^ Tr V BKP'BKP^ 

y \^\ ji+j2+2=j 
1 J 



E YAG,v,w). (14) 

l*-^l v,weG 

Note that the v = w terms vanish by ([H]). Each Yj{G,v,w) with v ^ w has zero mean 
because different rows of B are independent, with entries having zero mean. Thus, T{z^) 
has also zero mean. Again because of independence of rows of B, when we compute second 
moments, the terms in the sum below with {v, w} ^ {v' , w'} vanish. That is 

E{Tg{z')Tg{z^)) = ^ \2 ^Y,{G,v,w)Y,{G,v',w') 

v,w,v ,w SG 

= ij^ E ^Y,iG,v,w)[Y^iG,v,w) + Y^iG,w,v)] 

= I E ^Y,iG,o,v)Y^iG,o,v) 
veG\{o} 

For the last equality we fixed a vertex o of G and used the vertex-transitivity of the graph 
and the symmetry of Y in its last two parameters. ■ 



Lemma 12. For G finite, \TG{f)\ < cl\G\^l^ 

Proof. This follows directly from ©, dUD, and ([I2D. ■ 

We are now ready to prove Theorem [TJ 

Proof of TheoremUl Lemma[7]and Lemma [T2] show that E[Tg(-)Tg(-)] and Hg{-, ■) are || ■ ||*- 
continuous in /, g. We conclude the proof by approximating / and g with polynomials and 
using Lemma E] and Lemma [TTl ■ 
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5 Asymptotic normality 

For this section we consider a sequence of transitive graphs G„ converging to locally to limit 
G. The assumptions from the introduction hold for all G„ and the limit G with the same 
constant Ci. We may also assume that all graphs have degree d. 

We will prove convergence of certain sequences to normal under a topology which we 
introduce now. 

Let y denote the space of probability meas ures on R with finite second moment, equipped 



with the 2-Wasserstein distance, call it d2 (see IVillanil . Chapter 6). For two measures fi, v in 



the space, this is the minimal L^-distance over all possible couplings of them, i.e. 

1/2 



.(..:^,.(/,.-.P ..(.., 



where the infimum is over the set of probability measures on M^ with first marginal /i, and 
second v. 

The space (3^, d-i) is complete, and its topology is stronger than weak convergence. It is 
easy to show that for a sequence X„ of random variables, X„ — > X in this topology (i.e., the 
corresponding laws converge) if and only if X„ -^ X weakly and EX^ — * EX^. In particular, 
in this topology the function variance, Var : 3^ — > M, is continuous. 

We will need the following case of Lemma 2.4 of Chen and Shao (2004). It is a normal 
approximation theorem for dependent variables. It is proved with the use of Stein's method. 

Lemma 13 (Chen and Shao (2004)). Let I he an index set, {Xi -.i El} a family of random 
variables, and for A dX, let X^ '■= {Xi : i E A}. Assume that 

(1) For each i El, there exist Ai G Bi G Gi G X such that Xi is independent of Xx\Ai , 
XAi is independent of Xj\Bi , and Xb^ is independent of Xx\Ci ■ 

(2) There exists a constant 7 so that for all i eX we have 

max(|iV(Q)|,|{j:zeC,}|)<7, 

where 

N{Gi):={3EX:Gir\B,^^. 

(3) Each Xi has zero mean and finite variance, and W := ^jgjXj satisfies YaxiW) = 1. 
Then for 2 < p < 3, we have 



11 



sup\F{z)-^{z)\ <757^-^ VE|Xi|P, 

where F, $ are the distribution functions of W and of the standard normal N{0, 1). 

An immediate consequence is the following convergence result. 

Lemma 14. For any real polynomial f , the sequence T^Xf) converges as n —* oo in the 
2-Wasserstein distance to a normal random variable with zero mean and variance Hdf, /)• 

Proof. We will apply Lemma [T3l Let f{x) = X]j=o'^i'^"'' ^^^ -^ ^^^ ^^^ °^ vertices of G„. In 
the following, we will omit the subscripts for the matrices M, B associated with the graph 



Gn- We have 



\y^n\ \j={\ ki+k2+2=j J veGn 

where 

k . 

V\Gn\ j=o ki+k2+2=i 

It holds lim„_^_|_oo Var TcXf) = Hci^f^f) by Lemma [TOl If this limit is zero, then the 
sequence Tc^if) will converge to 60 in the 2-Wasserstein topology, and the result is proved. 
We may therefore assume that the limit is positive, and Var Tc„{f) > for all n. 

For f G X, define 

""^ - Var TaMY ~ k. " Var TaMY 

and the sets 

Ay : = {w : dist(f,w) < A;}, 
By : = {w : dist(f,w) < 2A;}, 
Cy : = {w : dist(f , w) < 3/c}. 

These sets satisfy the conditions of Lemma [13] (because the X„'s corresponding to vertices 
that are distance at least fc+l apart are independent random variables), and for all f , we have 
|A^(C„)| < |{w:dist(t;,w) <hk}\ <d{d-lY''-\ and \{w: v e C^}\ = \Cy\ < d{d - lY''-\ 
[Here we used the regularity of the graphs, but of course a uniform bound on the degree of 
the vertices of all graphs would work in the same way.] Also, the X„'s have zero mean and 

12 



finite variance by Lemma [TTl Pick any p G (2,3]. Lemma [T3l gives that the distribution 
functions F^r, Fz of W and of the standard A^(0, 1) satisfy 

sup \Fw{x) - Fz{x)\ < C\Gn\^~^E\A\P (15) 

xeR 

with C := 75{d{d - i)5'^-1)p-i and 

k . 
i=0 ki+k2+2=j 

which is independent of n for large n because the sequence converges locally to G. Thus 
relation flT^ . IG^I -^ oo, and p > 2 imply that the sequence 

Var TgM) 

converges to a standard normal random variable. By Theorem [1] and Lemma [10] we have 
Var TcXf) — > Hcif, /), and the result follows. ■ 

Our next goal is to strengthen Lemma [T^ to get Theorem[2]by showing that its conclusion 
is true also for all functions f in X. The proof is simply by approximation. We will use the 
fact that the set of polynomials is || • |L— dense in X (Lemma E]) and the following simple 



lemma (Lemma 38, Chapter 7 of iRoydenI (119881 )). 



Lemma 15. Let X be a metric space, and Y a complete m,etric space. Assum,e that fn '■ 
X -^ Y is an equicontinuous sequence of functions that converge pointwise on a dense subset 
of X . Then the sequence fn converges pointwise on the entire X to a continuous limit. 

We are ready to prove Theorem [21 

Proof of Theorem\^ Let T„(/) denote the distribution of Tg„(/). The sequence of functions 
Tn : X ^ y is uniformly equicontinuous because, by Lemma [7] and Theorem [T] we have 

nTGAf)-TGM\'<4\\f-9\\l 

By Lemma [HJ they converge pointwise at polynomials, which form a dense subset of X by 
Lemma [6] By Lemma [TSj, the limit T(/) of Tn{f) exists for all functions f E X and is 
continuous. Also, for the 2-Wasserstein distance between T(/) and T{g), we have 

d2{f{f),f{g))<cl\\f-g\U. (16) 
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Since the limit is normal on a dense set of points, and limits of normal random variables are 
normal, it follows that all limits T(/) are normal. Also the functionals 

/^Varf(/), f^HcifJ) 

are || • ||*- continuous (the first because of ( IT6l) and a property of the 2-Wasserstein distance, 
the second by Lemma[7]) and they agree on a || ■ ||*-dense set by Lemma [T^ Thus, VarT(/) = 

HG{f,f)ioTa\\fex. m 

6 A formula for the covariance 

The goal of this section is to reduce the problem of computing the covariance kernel (3 of 
Theorem [3] to inverting a Stieltjes transform. 

When the limit graph G is bipartite (see the discussion before Theorem [3]), we would like 
to write the limiting covariance in the form 

HgU os,gos)= [ f'{x)g'{y)d(3{x, y), (17) 



where s{x) = x"^. It is sufficient to check the identity flTTl) for the functions / = px^ and 
g = j9^2, since their linear span is a dense subset in X with respect to the norm || ■ ||^, (Lemma 
E]), and the integral on the right is clearly continuous with respect to the product topology 
based on this norm. For this choice of /, g the right hand side of 0171) equals 

/ (1 -x^xy (1 -^^yy ^^^"-^^ = l^^^^'-I ix-\-^) (y-V^) ^^^"'^^- 

Note that p\—px20s is an odd function. For bipartite graphs G it is clear from the expressions 
f ITU]) and fITT]) that Hq vanishes if one its arguments is odd. So we have 

^g(Pa2 o s,p^2 os) = Hcipx^p^). 
In light of the expression (j3]) we arrive at the following: 
Proposition 16. // the Stieltjes transform relation 



/"I 1 
7 T^T ^dp{x,y) {U 
„^„ i^ - ^ ) (y - f^ ) 

holds for all |A|, |/i| < 1, then (IT7|l holds for all f,gEX. 
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7 Explicit formulas for trees 

In this section, we look at the case where the hmiting graph G is the infinite d-regular tree 
Trf, we compute the measure d(3 introduced in (fTTj) . and thus complete the proof of Theorem 

El 

From now on we define square roots for complex numbers with a branch cut at the 
negative real axis. More precisely, for z = |z|e*^, with 9 E {—n,n], we set ^/z = ^J\z\e^^l'^. 



Proof of the covariance formula in Theorem We start by computing the left hand side of 
(ITSl) explicitly. The Green's function for the infinite c?-regular tree, Q\{v, w) = ^A(dist(f , w)), 
is given by 

Gxir) = bxa\, 

where 

1 - V^l - pA2 ^^ 

"^^=^ 2(^-l)A ' ^^-(1-^^^) ' 
and p = 4{d — l)/(f' as in the introduction. See IWoesd (120001 ) . Lemma 1.24. 



First, an observation. For two vertices i,j, let i denote the unique neighbor of i closest 
to j, r the distance between i and j, and define B*- = Ba. For A with modulus less than 1, 
by symmetry and since J2k ^ik = 0, we get 

(BQ^),, = B^eGxir - 1) + ^^..^^(r + 1) = 5*,.^*(r), 

where 

Qlir) = Qxir - 1) - Qxir + 1) = 6a(1 - aD^"'- 
Thus for any vertex w ^ a, with r = dist(o, w), we have 

E{Bg,)^,UBg,)o,u. = EB:lgi{r)g;{r) = gi{r)g;{r), 

because E_B*^^ = 1, and therefore the sum in the left hand side of flTSj) equals 

E 2gi{rrg;irr = 2dJ2{d-ir-'giirrg;irr = 2 ^ i-(d-i)(a a / • ^^^) 



Now this can be expressed in terms of s = a/1 — pA^ and t = a/1 — pp^. Indeed, in these 
variables, we have the simpler expressions 



2 _ 1-s L _ 2(d- 1) 

(rf-l)(l + s)' ^~ {d-2) + sd' 



and (fT9|) becomes 

32{d-l)d 



(1 + s)(l + t) {{d - 2)(1 + s)(l + t) + 2(s + t)) ■ 



(20) 



Introduce new variables u,v by u = \ ^,f = /i ^, and let P{u,v) denote expression (120|1 as 
a function of them. In terms of these variables, relation fllSp becomes 



uv J {x — u) {y — v) 

Due to our convention for square roots (see beginning of this section), the quantity s is an 
analytic function of m in C \ [0, p], and the same holds for t as a function of v. 

Assume that d > 2. Then, the denominator in fl20l) does not vanish because s,t have 
nonnegative real parts (if we set the last factor in the denominator equal to zero and solve 
for t, we get a quantity with negative real part). So /3 defined by fl2Ul) is a holomorphic 
function of {u,v) on (C \ [0,p])^. In fact, even the limits of the denominator when m or f 
approaches [0,p) are not zero. 

Since the function h{u,v) := {uv)~^l3{u,v) is holomorphic in (C \ [0,p])^, and it decays 
as {uv)~^ near infinity, Cauchy's formula will express its values in terms of double contour 
integrals around the segment [0, p]. Shrinking the contour around [0, p], we get a line integral, 
and we take into account the different limits of h as one of its arguments approaches the 
segment from the upper or the lower half plane. That is, when u approaches x G [0, p) from 
the upper half plane, we have s{u) — > s{x). While when u approaches x from the lower half 
plane, we have s{u) -^ —s{x). The difference comes from the branch cut discontinuity of 
square root in the definition of s. Thus we have 

(5{x,y) I f fp 1 (3{u,v)dudv 



xy 47r^ Jo ^o {u — x){v — y) uv 

where for m, f G [0, p] we have 

(3{u,v) = y, OT(3[oS,Tt] 

(T,r=±l 

hl2{d-lfd'^{d-2)st 
{d\s^ + t2) -{d- 2)2(1 + sH^)f - (8(d - l)sty 

and [3[as,Tt\ refers to the expression (l20l) with (s,t) replaced by ((Ts,rt). 
Consequently, the density of the measure d/? is 

1 (3{u,v) 

f3d{u,v) = -—^ • 

4tt^ uv 
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Substituting the expressions of s, t in terms oi u,v, we find 

Pd{u, V) = -^d [d - 1) (d - 2) . 

with 
A : = mV ({d^is" + t^) - (rf - 2)2(1 + s¥))2 - (8(rf - l)sty 

= p2 16(2d - 3)(m - vf + {d- 2)2 /'pfi:(!i±^)2 + 4(rf + 3)(m - t;)^ + p2(rf _ 2)'' 

where k,{u) = 2d^Ju{p — u). Then 

^ , , 2/«(M)KM(rf-2)d4 7r-2 

PdiU,V) = T) . 

16(2rf - 3)(m - w)2 + (d - 2)2(pA;(^) + 4(rf + 3)(m - t;)2 + p2(rf _ 2)2) 

This proves the d > 2 case. 

The case d = 2 can be easily shown by using continuity in the formulas as d I 2. We get 

/32(a;, y) = - K,{x)6^{y) = —\/x{l - x) 6^{y). 

TV IT 

The qualitative difference here is that for d = 2 the denominator of ( l20|l does vanish along 
a line. ■ 

Acknowledgments. This research is supported by the Sloan and Connaught grants, the 
NSERC discovery grant program, and the Canada Research Chair program (Virag). We 
thank Amir Dembo for encouraging discussions. 

References 

G. W. Anderson and O. Zeitouni (2006). A CLT for a band matrix model. Probab. Theory 
Related Fields, 134(2) :283-338. 

A. S. Deif. Advanced matrix theory for scientists and engineers. Abacus Press, New York, 
second edition, 1991. 

P. Diaconis and S. N. Evans (2001). Linear functionals of eigenvalues of random matrices. 
Trans. Amer. Math. Soc, 353(7) :2615-2633. 

A. Gamburd, S. Hoory, M. Shahshahani, A. Shalev, and B. Virag. On the girth of random 
cayley graphs. To appear in Random Structures Algorithms. arxiv:0707.1833. 

17 



R. A. Horn and C R. Johnson. Topics in matrix analysis. Cambridge University Press, 
Cambridge, 1991. 

A. Klein (1998). Extended states in the Anderson model on the Bethe lattice. Adv. Math., 
133(1):163-184. 

N. Linial and A. Widgerson. Expander graphs and their applications. Lecture notes available 
at http://www.cs.princeton.edu/~boaz/ExpanderCourse/. 

B. Rider and B. Virag (2007). The noise in the circular law and the Gaussian free field. Int. 
Math. Res. Not. IMRN, (2): Art. ID rnm006, 33. 

H. L. Royden. Real analysis. Macmillan Publishing Company, New York, third edition, 1988. 

C. Villani. Optimal transport, old and new. Available at http://www.umpa.ens- 
lyon.fr/~cvillani/surveys.html. 

J. vom Scheldt and W. Purkert. Random eigenvalue problems, volume 63. Akademie-Verlag, 
Berlin, 1983. 

W. Woess. Random walks on infinite graphs and groups, volume 138 of Cambridge Tracts in 
Mathematics. Cambridge University Press, Cambridge, 2000. 

O. Zeitouni. Lecture notes on random walks in random environment. In Lectures on proba- 
bility theory and statistics, volume 1837 of Lecture Notes in Mathematics. Springer- Verlag, 
Berlin, 2004. Lectures from the 31st Summer School on Probability Theory held in Saint- 
Flour, July 8-25, 2001, Edited by Jean Picard. 

Dimitris Cheliotis. Eindhoven University of Technology, Eurandom. L.G. 1.26, P.O. Box 
513, Eindhoven, 5600 MB, The Netherlands. dimitriscOgmail.com 

Balint Virag. University of Toronto, 40 St George St. Toronto, ON, M5S 2E4, Canada. 
balint@math . toronto . edu 



